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Abstract 

General solution of RG equations in the framework of background perturbation theory is 
written down in the large N c limit. A simplified (model) approximation to the general solution 
is suggested which allows to write /3(a) and a(/3) to any loop order. The resulting oib{Q 2 ) 
coincides asymptotically at large \Q 2 \ with standard (free) a s , saturates at small Q 2 > and 
has poles at time-like Q 2 in agreement with analytic properties of physical amplitudes in the 
large N c limit. 
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1. The running a s of Standard Perturbation Theory (SPT) has unphysical singularities (Landau 
ghost pole and cuts) in the Euclidean region; which contradict analiticity of physical amplitudes 
and are not observed in lattice calculations of a s (Q 2 ) Q. 

As it was shown in [| the Background Perturbation Theory (BPT) is free of these defects 
and the coupling constant as of BPT displays an important property of IR freezing (saturation) 
with ccb(Q 2 — 0) ~ 0.5 ||, B|. This behaviour of is well confirmed by experimental data on 
spin splitting of quarkonia levels ^, [| and by lattice data || . 

Analytic properties of asiQ 2 ) in the Q 2 plane were not studied heretofore, and is the topic 
of the present letter. To this end the large N c limit is exploited which ensures that any physical 
amplitude has only simple poles @|. This requires that as(Q 2 ) should be a meromorphic function 
of ()'. 

To find it explicitly we study below the general solution of RG equations and suggest a simple 
model approximation to it, yielding /3(a) to all orders and explicit form of as(Q 2 ). 

Resulting as(Q 2 ) can be compared with a s (Q 2 ) of SPT and coincides with the latter asymp- 
totically at large Q 2 . 

2. Separating background field in the total gluonic field = B^ + a M and expanding 
in ga^ one gets the BPT series, where coefficients depend both on external momenta and the 
background field gB^, which turns out to be RG invariant, if background gauge is chosen ||. 
Correspondingly as = <? 2 /47r satisfies the standard RG equations, where momenta and averaged 
background characteristic tr < F fJil „ 1 {B)F^ 2 v 2 (B)...F^ „ (B) > enter together and will be denoted 



One starts with the standard definition 

din /i 



PM (i) 



and represents /3(q.b) in the form 
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where if(z) is an arbitrary function with some conditions to be imposed below, and prime means 
derivative. 

Solving ([!]) with the help of (||) one has 

47T 

aB = /3[lnM 2 +X+|M^)] 

Here x is an arbitrary dimensionless function of {V 2 }, while <p is partly fixed by few coefficients 
of expansion of /3(a), 

q{ \ A 2 A3 ^2 4 lA s 

3(a) = -—a - — —a° - ——^a - (4) 
y ' 2n An 2 6Att 3 (47r) 4 v ' 

(Note that [3i are defined as in for nt = 0, N c = 3 and MS scheme 0o = 11, /3i = 51, P2 — 
2857, (3 3 = 58486.) 

At this point we shall use the large N c limit and require that as(Q 2 ) be a meromorphic function 
of Q 2 . Consider e.g. the process e + e~ into hadrons, and the photon self-energy part n(Q 2 ), which 
has a pole expansion ||] 

U(Q 2 ,a B ) = £ n2 C ; l{ ^ , + subtractions (5) 
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Expanding in as and doing a partial summation (equivalent to the summation of leading 
logarithms) one obtains a running as(Q 2 ) which has poles at some generally speaking, shifted 
positions as compared to the poles of the lowest order - M 2 (0). The latter are proportional 
to n; M 2 (0) = rain + M§; n = 0,1,2,..., m 2 , = Atto, C„(0) = m 2 , |[. The lowest (partonic) 

approximation thus reduces to n(Q 2 ,0) ~ ip( — ), where ip(z) = thj ■ 

To specify the form of the function and x(Q 2 )i we shall use two requirements: 

i) the correspondence principle with the SPT, which tells that for Q 2 3> k 2 , where k 2 is the scale 

of nonperturbative vacuum fields, our solution should coincide with that of SPT. 

ii) ccb(Q 2 ) should be meromorphic function of Q 2 , and hence both x an d should be meromorphic 

functions of their arguments (Q 2 and respectively). 

The first requirement means that coefficients C n in the sum over poles should have finite limit 
Coo, i-e. one can represent this sum as the Euler function ip(z) — ^r§\ plus a sum with fast 
decreasing coefficients, which we write symbolically as a finite sum. Thus the general solution 
satisfying i) and ii) can be written as 

hl/i +X = ln A 2 - + ^ ( ^ 2 — ^T.qT^ (6) 

n=l " 

¥>(— )=^(— + A)+pf— ) (7) 
olb a B \ttfl / 

where A is some constant and A = Att/Pq. Note that our solution (||), (|7|) satisfies (jl]) and 
reproduces the first two coefficients in the expansion (Q) when all b n = and P = 0, while our 
02, 03 appear to be smaller; however the latter are scheme-dependent. It is easy to reproduce p2, 
/?3 by choosing in (Q) P ^ 0. This choice of <p(—) leads to a small numerical (around 2%) change 
in as as compared to the minimal choice (||) for all Q 2 < 0. 

3. In what follows we confine ourselves to the minimal model of (||), (R) with all b n , and P 
identically zero, which yields 



4-7T 

Po 



'. m 2 , ,Q 2 + M 2 2/3! , . A 



(8) 



Eq.(||) defines as as a meromorphic function of Q 2 in the whole Q 2 complex plane. 

Moreover, Eq.(||) defines as to all orders in the "loop expansion" obtained by iteration of 
the last term in the denominator of (||) and will be compared below to the SPT loop expansion. 

Consider e.g. large Q 2 , v "T " > 1- Then using the asymptotic expansion for z — > oo 

,(,)_,„-£-£;*», (9) 
fe=l 

where P„ are Bernoulli numbers, one obtains to the lowest order 

a * " J* < 10 > 

which explicitly shows the AF behaviour for large Q 2 and absence of Landau pole (for Mq > A 2 ). 

The forms (^) and (||) admit in general poles of as(Q 2 ) for Q 2 > and poles of /3(ob) for 
as > 0, which would be unphysical. The latter poles are defined by equation 1 = J^Bo ^' (a + A), 
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and are excluded from the region a > when A > Ao = 1.255. This latter condition also excludes 
singularities of ub{Q 2 ) (@)j at Q 2 > 0. At the same time /3(a) has a point of condensing zeros on 
the negative side of a = 0, and /3(a) is monotonically decreasing when a — » oo, as it is shown in 
Table 1. Note that a pole of /3(a) at a = a* > leads to a situation when real solutions a(q) exist 
only for a < a* , Q > Q* . 




Figure 1: Exact solution of Eq.(8) with A = 1.5 as a function of s = Q 2 (solid line), and exact 
solution for as from the extended model, Eq.(3) with tp from (21) and rj = 3, A = 1.5 which 
reproduces all known (3 n coefficients (broken line). 



4. Although chb(Q 2 ) does not have poles for A > Ao in the Euclidean region, Q 2 > 0, it has 
an infinite number of poles for Q 2 < 0. 

These poles are given by zeros of the denominator in (§1), which can be written as 

- Q 2 = s n = a® - S n , 4 0) = M 2 = Ml + nm 2 (11) 

and S is then found from (||) to be |l^| 

M 2 

S-'-ln^. (12) 

Hence poles of as are shifted as compared with poles of ip ( m2 ) , the latter describe the 
physical states of (double) hybrids - the background analog of gluon loops renormalizing a s in 
SPT. In the vicinity of the pole as has the form |[o) 

(47T/A, - V/(A)) M, 2 - s ( 

Hence perturbative series of BPT is a sum over poles of olb and nonperturbative poles of 
n(Q 2 ,0) in (||) (dual to partonic contribution), all poles being in the region s > Mq — 5 > and 
this situation is in accord with the large N c analytic properties. Now one has to show that the 
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"physically averaged value" of as(s) corresponds to the AF behaviour, as it is seen e.g. in the 
e + e~ annihilation. To this end one one can choose two equivalent procedures: i) to introduce the 
width of the pole in as(s), or ii) equivalently to consider the BPT series at the shifted value of 
s,s= s(l + ry) (see e.g.p|). 

In what follows we shall exploit the second path. Considering Imas above the real axis, at 
s(l + £7) with s large and 7 fixed 7 <C m allows to use asymptotics ( |l0|) (in one loop order) with 
the result 

■^"'^"■ ^igr^-^ (i4 » 

Another form of discontinuity results from the expansion of the Adler function D(Q 2 ) = 
Q 2 d J\^)i , from which the hadronic ratio R(s) can be obtained as |l2| . 

-s+ie 



R(s) = ± r* e D(*)^ = ± [ D(a)^ (15) 
2m J_ s _ ie cr 2m J c{s) a 



Here C(s) is the circle of radius s around the origin, comprising all singularities of D(a) in 
Euclidean region. 

Similarly to [fl2| one can define the operation transforming the known perturbative 

series for R(s); 



R(s) = N c J2e 2 q 
q 



1 -!>{ — } • 



where 



V 7T 



fOB_\ k \ _J_ f ( Qa(g) \ fc dcr 



(16) 



We shall consider the situation, when one can use the logarithmic asymptotics ([!(]) for olb 
in the Minkowskian region. This happens in the physical situation described above: when poles 
have imaginary parts (widths of resonances), or equivalently when one consider effective R(s) = 

5j(II(s(l + £7)) — II(s(l — 17)) with fixed 7 (so that requirement |arg(z = M ^ 2 S ) < 7T is fulfilled). 

In this case one obtains, inserting in the integral ( |l7| ) the asymptotic form ( |l0| ) (sec [[Io| for 
details) 

f «s 1 4 f 1 7T 2M 2 / M? \ 1 

* i — > = t- \ -arctg- -pr + — + O 2-2 [ 18 

I 7T / 6 W S lns/A 2 S ln 3 s/A 2 \sln 4 s/A 2 J} V ; 

5. One can now compare the power ("loop expansion") of as in (|^) in powers of ar2 ( |Io|) with 
the SPT expression for a s [0, O. To this end we shall write it as follows: 



(0) J 2A lnL 4/3 2 
where L = ln( A2 ) and for n/ = 0, A = 1.5 



dni-i) 2 



(19) 



& = -|L(A-i)-i = -1.436 (20) 

This should be compared with the MS value of SPT, bjjg = 0.26 (rif =0). 

To bring our theoretical expressions for /32, /?3 in agreement with the computed MS values (see 
|| for the corresponding references), one should generalize the minimal model considered above 
by adding a term to ^(~): 

^(I) = AV/(- + A) + -^ (21) 
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Choosing C = 1.253 and arbitrary real 17, one obtains both /3a, [}% in close agreement (within 1%) 
with MS values. In this case also b in ( pp| ) gets contribution from the new term and becomes 
b = b + '^jj^-C = 0.263, which agrees with the SPT ||. We study numerically both minimal 
model and the corrected one as explained above and present the results for /3(a) in Table 1 and 
for a B {Q 2 ) in Table 2. 

From Table 1 it is clear that nonperturbative /3(a), namely (3 m i n {a), calculated as in (0), (0) 
with P = and j3 ex t{a) calculated as in (ph, (f2l|), both decrease more moderately (as (-const a 2 )) 
at a > 1 as compared to purely perturbative 3 loop expression (Q). This might indicate that the 
series (|4|) is asymptotic and should be cut off after few first terms (note that P% is larger than 
previous three coefficients). 

For Q 2 = (10 GeV) 2 and Mq — 1 GeV the correction in as(Q 2 ) introduced by replacements 
( pl| ) is around 2% and it decreases for Q 2 growing. 

Defining A = A# = Aj^s (since olb and a s practically coincide for Q 2 > 40 (GeV) 2 ) one 
can compare exact as(Q 2 ) obtained by solving ^ with a s (Q 2 ) (two loop) of SPT and observe a 
significant difference only at small Q 2 , Q 2 < 10 GeV 2 . 

One can see in Table 2 that the difference of a^™ m ' ) and a~ 2loop ^ is 20% at 1 GeV 2 and drops to 

10% at 5 GeV 2 , which implies that the use if a^™ m ' ) instead of a s in the Euclidean region would 
not produce much change in the region Q 2 > 1 GeV 2 while it improves the situation for Q 2 < 1 
GeV 2 , where a^" saturates in agreement with experiment and lattice data (see fl|]-|§ and (m) 
for the discussion of this point). 




Figure 2: f3 m i n (a) in the minimal model Eqs.(2), (7) with P = 0, A = 1.5, reproducing scheme- 
independent coefficients /3o,/3i- dashed line, and /3 ea; t(a) in the extended model Eqs.(2), (21) with 
e = 3 and A = 1.5 (solid line) and A = 3 (dotted line). 

6. We have used the large N c limit to simplify analysis of analytic properties of background 
coupling constant as- This has enabled us to study as(Q 2 ) also for time-like Q 2 and to suggest a 
simplified model solution of exact RG equations, which yields ag to all orders and prescribes the 
nonperturbative expression for /3(a). 

The resulting solution for coincides with the standard a s at large \Q 2 \ and displays fa- 
miliar AF properties. At the same time ub{Q 2 ) stays finite everywhere in the Euclidean region, 



(i 



demonstrating phenomenon of IR saturation. For time-like Q 2 the model as(Q 2 ) in the large N c 
limit acquires pole singularities, which have been treated above in the averaging procedure. It is 
worth noting that our minimal model correctly reproduces both scheme-independent coefficients 
0o,0i- At the same time scheme- dependent coefficients 02,03 etc. are defined by additional pole 
terms with arbitrary constants in them. Several improvements are possible with the suggested 
model solution. First, one can consider rif ^ 0, N c finite, in which case poles would acquire 
width. This will not change much the averaged values of a_e(s) calculated above, but will enable 
one to make local predictions of resonance behaviour. Secondly, our model is too simplified, since 
ip( m 'i ) appearing in fl8J), is responsible only for the equidistant spectrum of hybrids, (equivalent 
of gluon loops), but this spectrum should not be equidistant for lowest states, which yields power 
corrections, considered in p~C| ] . 

The author is grateful to N.O.Agasian and A.V.Radyushkin for helpful discussions; a partial 
financial support of RFFI grants, 00-02-17836, 00-15-96786 and INTAS grants 00-00110 and 00- 
00366 are gratefully acknowledged. 



a 





0.1 


0.2 


0.5 


1 


1.2 


1.5 


2 







-0.0187 


-0.08 


-0.587 


-2.85 


-4.36 


-7.38 


-14.59 


0ext{a) 





-0.019 


-0.082 


-0.65 


-3.36 


-5.23 


-9.04 


-18.43 


f MS 





-0.01895 


-0.0833 


-0.767 


-6.82 


-13.56 


-33.35 


-115.24 



Table 1: The QCD /3-function 0(a) for nf = 0, iV c = 3 obtained numerically to all orders in a for 
the minimal model m in(a) Eq.(2) with tp = ^ + A, and A — 1.5; A ~ 1.14, and for the extended 
model ex t with if' from (21). The bottom line refers to 0(a) in MS as given in Eq. ([I]). 



Q 2 (GeV 2 ) 





0.5 


1 


1.5 


2 


3 


5 


10 


40 


100 




0.546 


0.394 


0.337 


0.306 


0.286 


0.26 


0.232 


0.201 


0.159 


0.139 


(ext) 
«B 


0.754 


0.462 


0.379 


0.337 


0.311 


0.278 


0.245 


0.210 


0.163 


0.142 


(2loop) 

a B 


0.64 


0.359 


0.328 


0.293 


0.27 


0.244 


0.217 


0.188 


0.15 


0.131 


(2loop) 

a\ 




0.73 


0.404 


0.33 


0.293 


0.256 


0.229 


0.19 


0.149 


0.131 



Table 2: as(Q 2 ) computed numerically to all orders within minimal model, Eq.(||), with A = 0.37 
GeV, A = 1.5, m = M Q = 1 GeV and in extended model, Eqs. (§), @, with rj = 3, C = 1.253 
reproducing all known coefficients /3„, n = 0, ...3 in MS scheme for nj = 0. The last two lines refer 
to the two-loop approximation of as and a s (in SPT) with the same parameters. 
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